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DISCLAIMER 



This report was prepared by consultants under contract 
to the Air Resources Branch, Ontario Ministry of the Environ- 
ment. The views and findings expressed herein and the quality 
and accuracy of the text are the sole responsibilities of the 
contractor, Opinions and recommendations expressed are those of 
the contractor and should not be construed to represent policy 
of the Ministry of the Environment or the Government of Ontario, 
Mention of specific brand or trade names does not constitute an 
endorsement by the Government of Ontario. 
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SUMMAKT 

A mathematical model is derived which enables evaporation 
rates to be calculated for spills of liquids, such as PCBs, of low 
to moderate volatility, and permits estimation of atmospheric 
concentrations of the volatilizing substance at points in the 
vicinity of the spill site. The model is fitted to data obtained 
from the vicinity of a site contaminated by PCBs . The model is also 
applied for illustrative purposes to an actual PCB spill, A simple 
calculation procedure is defined "which can be used for on-site 
estimation of exposure hazard to evaporating substances, 
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1, INTRODUCTION 

When a chemical substance is spilled on land, evaporation may be 
a significant mechanism of mass loss. A knowledge of the evaporation 
rate is useful in several respects. For volatile liquids (such 
as gasoline) it will determine the "lifetime" of the spill. For 
flammable substances it will influence the extent of fire hazard 
in the vicinity. For toxic substances it will influence the 
dimensions of the zone in which it is necessary to take protective 
measures to avoid excessive exposure . 

It is thus valuable to have the capability of calculating, 
at the spill site, the rate of evaporation and the likely atmospheric 
concentrations at various points in the vicinity, A desirable 
feature of such a capability is that the calculation can be done on-site, 
by the personnel involved, without the need to use sophisticated computer 
programs or have access to advice from others. In this work we devise 
and use such a calculation procedure in the full expectation that it 
will be subject to future improvement, The primary aim is to 
establish the framework for a reasonably simple procedure which is 
based on physically reasonable assumptions . 

It is assumed that the person undertaking the calculation has 
available the following information: 

(i) chemical identity of the spilled substance, 

(ii) the substance's vapour pressure or at least its boiling point. 

(iii) the dimensions of the spill ( area, volume and thickness or 
depth) 

(iv) ambient meteorological conditions of temperature, wind speed 
and direction, 

(v) knowledge of the terrain topography. 



- 2 - 



For interpreting the results of the calculation it may be useful 
also to have information on acceptable exposure concentrations, eg. TLVs 
or flammability limits. 

The approach taken is to devise in Section 2 a mathematical model 
of the evaporation and dispersion process. In Section 3 this model 
is fitted to experimental data obtained from a site in Ontario in which 
there had been a spillage of poly chlorinated biphenyls (PCBs) . In Section 4 
the model is applied to a second case in which a quantity of PCBs was 
spilled in a school yard as a result of transformer failure and in which 
air concentration data were also obtained. Some conclusions and recommendations 
and discussion of limitations are given in Section 5. A, draft "manual" 
is included in Appendix I for undertaking the calculation. 
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2. MATHEMATICAL MODEL 

2.1 Approach and Limitations 

In this section we first address the question of calculating the 
evaporation rate of a spill. This is followed by derivation of a plume 
dispersion model. The two are then linked. This linkage is simple when 
the substance is fairly invoiatile and the concentrations achieved 
over the pool are well below the saturation values, The rate of 
evaporation is then little affected by the concentration in the overlying 
air and the evaporation rate is kinetically controlled. In cases when the 
plume is advected only slowly from the pool there may be "blaaketdtng" 
and the atmospheric concentration reduces the evaporation rate. 
Finally there is a brief review of vapour pressure determinations for 
liquids in general but for PCBs in particular. 

It is assumed that the pool temperature is known. This will 
generally be close to ambient air temperature but there may be evaporative 
cooling and heating by solar radiation; thus, some judgement is necessary 
to include these effects. 

If the liquid is absorbed into a porous substrate such as soil, 
the absorbed liquid is shielded from evaporation and evaporation will 
be retarded. This effect is not included in this work; although 
methods of addressing the issue are discussed later. 

Throughout this work SI units are used, but conversion factors 
are given in Appendix I , 

2.2 Evaporation Rate 

When a quantity of liquid is spilled on the ground, it rapidly 
achieves environmental temperatures and becomes subject to spreading, 
evaporation and penetration into the underlying ground. If the 
surface is water-logged, there will be relatively little penetration 
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but spreading will be enhanced . 

If a volume Vm^'ia spilled over an area Am 2 , the depth is D or V/A m. 
It will evaporate at a rate controlled by the liquid's vapour pressure 
P (Pa) (as influenced by temperature) and composition (mole fraction x) 

and the evaporation mass transfer coefficient K , which generally has 

-2 
a value of appoximately 10 m/s at normal wind speeds (Mackay and 

Matsugu, 1973) . Specifically the rate of evaporation Q g/s will be 

Q = K APxM/(RT) 

e 

3 
where R. is the gas content (8,314 Pam /mol K) , T is absolute temperature 

(K) and M is the molecular weight (g/mol) * 

A convenient expression for K is a simplification of that devised 

by Mackay and Matsugu (1973) ; namely, 

n • 7ft 
Ke - 0-0025 U 

where U is the wind speed (m/s) « 

In many cases the liquid is pure; thus, the mole fraction term 
x is unity and can be ignored. 

In an askarel (PGB transformer oil), x is in the range 0*5 to 0*75, 
In practice, it is acceptable to interpret x as the volume fraction of 
the substance of interest. 

The equation for Q is useful in two respects. First, it gives the 
source strength of the spill which is used later in the dispersion 
calculation. Second, it can be used to estimate the spill lifetime t . 
Evaporation will be complete when Qt equals the mass in the spill (g) , 

ie. 10 Vd where d is the liquid density (g/cm ) , there being 10 cm 

3 
in a s , It follows that 

6 
T m 10 Vd/Q seconds 

= 10 6 Vd/[KAPxM/(RT)] 
e 

= 10 6 Dd RT/(K PxM) 
e 
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where D is depth, ie. V/A. 

Typically d is 1-0, T is 288K, K is 0,01 m/s, x is 1.0 and M is 



200 g/mol; thus, 

9 
T ■ 1.2x10 D/P seconds 

- 3^3 x 10 D/P hours. 

-3 
For a 1 mm. thick spill (D=10 m) of a volatile liquid of vapour 

pressure 10 Pa ie. 0*01 atm or 7*6 mm Hg, T is thus 0*33 hours or 

20 minutes . 
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2.3 Florae Dispersion 

Assuming that a constant source strength or a concentration 
immediately above the evaporating pool is known and is independent of 
the atmospheric dispersion processes, it is possible to calculate 
the concentration at points downwind if certain assumptions are made 
concerning 

(1) steady state vs. unsteady state (time variance) 

(2) wind speed and direction constancy with time and position 
and as influenced by pool and adjacent terrain roughness 

(3) diffusivity constancy with time and position in three 
dimensions 

(4) interaction of the substance with the ground (ie, reflection 

or absorption) 

(5) buoyancy effects of the substance at high concentration 

In general, different models arise from different assumptions being 
applied, the more severe the assumptions, the simpler the model, 

In this work, we reviewed dispersion models, particularly the 
conventional models reviewed by Turner (1970) and the vapour cloud dispersion 
model described by Reijnhart et al. (1980), The latter model is very 
effective but unfortunately does not lead to a simple expression for 
concentration or mass flux. Accordingly, we have derived a novel 
model which is essentially an adaptation of the Gaussian model, which 
it is hoped will adequately describe the concentration profiles. It 
is derived here from first principles. 

If a point source is emitting substance at Q g/s into a wind 

of velocity U m/s and the plume spreads horizontally to width H and 

3 
vertically to width h, there being homogeneous concentrations (C g/m ) 

within these dimensions it follows that the conservation of mass 
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requires that 

Q - *CUHh (1) 

thus C = Q/ UHh (2) 

In practice C is not constant but varies in a near-Gaussian manner 
in both dimensions. Postulating that 

C = C exp(-y 2 /2g 2 - z 2 /2a 2 ) (3) 

y z 

where C is the centreline concentration gives the concentration at 

© - 

any width y or height z. Conservation of mass requires that the 
integrated concentration equals Q, y being integrated from - m to + * 
and z from Q to + <*» (if there is total reflection from the ground). Now 

</exp(-z 2 /2a 2 ) = a /Tv/2 (5) 

Q - C PCS /2tF)(c /2tF/2) (6) 

o v y z 

2 2 2 2 
oi C = (Q/Dtto a )exp(-y"/2a - z /2a ) (7) 

y z y z 



which is the conventional ground level emission equation. 

We are concerned with concentrations around the pool; thus we 
redefine y as the perimeter distance rather than the normal distance 
as shown in Figure 2.1. It follows that yean be derived now from 
an angle and the radius rather than a distance. The total perimeter 
at a radius X is 2ttX and is equivalent to 360° ■ thus, if 9 
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is the angle between the point of interest and the wind direction 

y = 92TTX/360 
Alternatively if 6 is expressed in radians § 

y m 0X 
We also redefine ay in terms of an angle a n as 

ay « O.2TTX/360 or <$>X 

Inserting these into the Gaussian equation gives 

C = Q(ttUg (a Q 2TTX/360))exp(-e 2 /2c Q 2 - z 2 /2a 2 ) 

z y w z 

or C - 180Q/(iT 2 Ua z a e X)exp(«9 2 /2a e 2 -z 2 /2a z 2 ) (8) 



If a becomes large there will be some 'ircraparQund" , ie. the 
substance will diffuse over 180° and will "double up" the concentration, 
In practice this does not occur appreciably if C7- is less than 60°. It 

y 
can be included by replacing 

exp(-e 2 /2a Q 2 ) by exp(-6 2 /2a Q 2 ) + exp(-(9 - 360) 2 /2a e 2 

This is similar in principle to the hypothetical "second stack" 
used to predict reflected concentrations. Effectively it doubles 
concentrations at ISO* 5 . In principle a (0 - 720) term could also be 
included but this is necessary only if a n exceeds 180°. In practice 
(as is shown later) if a n is 180° the concentration varies only by 
some 10% around the perimeter. There is thus no merit in using 
values greater than 180°. 

These concentrations are illustrated in Figure 2 .1. It is 
concluded that the difference in concentration around the perimeter is 
adequately expressed by the equation below provided that only values 
of a m less than 180° are used. 
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+m 
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-366 



- ICO -*• 



*360 
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As is shown later, when 8 is 180 the concentration at 180° is 
actually slightly greater than the concentration at 0°, which is, of 
course, physically unrealistic. This slight discrepancy can be readily 
avoided by setting a maximum value of 160° for 9 at which angle the radial 
distribution is quite even. 

C m 180Q/(Aa z a e X)exp(- z 2 /2g z 2 )[exp(-e 2 /2a Q 2 ) + exp (-(8-360) 2 /2o Q 2 ] 

If desired, the angular dispersion quantities cr fl and 8 can be 

replaced by the perimeter distances a and y, This may be advantageous 

j 

because it avoids any confusion between degrees and radians during 
calculation. Unfortunately, it is convenient to measure wind direction 
relative to the point: of interest as an angle. Considering these 
issues has led to a compromise in which <J is calculated and included 
in the equation, but a Q is also calculated from it. The calculation 
sequence is given later in the Appendix* 
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"Ehe centreline ground level concentration at 8 - 0, Z ■ 
becomes from Equation 9 

C - Q/CttuVj a )exp(- 2 2 /2ej 2 ) [exp(-8 2 /2g 2 ) + exp (-(8-360) 2 /2a 2 ) 1 

m y z 9 Q 

The concentration (Q/ttUo a ) is approximately the centreline ground level 

z y ° 

concentration. The actual value is slightly higher because of the term including 

including (8 - 360) which is not zero when 8 is zero. When or is 130°, 

6 

i.e. the "worst case", the difference is only 14%. 
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lt can be shown that this equation reduces essentially to 
Equation 2 under certain conditions. If a non-Gaussian expression is used 
for the vertical concentration, the concentration being assumed to be 
constant to height h t then 



or 



Q - C Uo /Mi 

f ° y 



C ■ Q/yziTlJa h i.e. h replaces a /S7/2 

o H -y z 



which leads to 



C - 180Q/ir 2 Ua XV^7ih[exp(-9 2 /2a^ 2 ) + exp(-(e-360) 2 /2a 2 ] 



When a» is 180 ° f the exponential terms on the right assume the 
following values for selected values of 8: 
8 - 0* , 1 + 0.135 * 1.135 
6 - 45* , .969 + 0.216 * 1*185 
■882 + -325 - 1.207 
755 + -458 * 1.213 
•607 + .607 + 1.214 
,191 and the maximum variation is only 0-056 



6 - 90° , 

9 - 135 , 
I - 180 f 
The mean value is 1 



or -4. 7%, 

The radial distance X is related to the perimeter H as H • 2ttx . 

Replacing X with (H/2tt) t the exponential group by 1.191 and G~ by 

y 

180° gives 

C = 0.95Q/UHh 
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which is essentially equation (2) derived earlier. It follows that 
Equation (9) can be expected to give concentrations as a function of 
radial distance X and height z which are consistent with the mass balance. 
The equation has been shown to reduce to the intuitively obvious 
Equation 2. . , „_ 

Correlations for q . q. and ay. 

The simplest procedure is to assume correlations equations in which 
q and q increase linearly with distance from the source. Specifically 

m Jr 

it is proposed that 

q * C-H + C, X 
z 3 1 

where H is the pool diameter, X is the distance from the pool edge and 

C, and C. are constants. Intuitively, C, should be small, reflecting 

3 1 j 

the height of plume rise over the pool. When X is zero,q becomes C-H. 



The dispersion coefficient diagrams for q under stability condition Q 



suggest that C^ is of the order of 0.07 m~ since a rises from 7.5 to 70 a 



z 

m e, 

K 

when X increases from 100 to 1000 m. C. will also be affected by terrain 

x 

roughness, trees, buildings, etc. which induce vertical flow. To a first 
approximation q will be increased by the height of the "object". 
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It is proposed that 

a - C,H + BJt/tt 

y 4 2 

where U is wind speed; thus, X/U is the diffusion time and C„ can be 

viewed as a velocity of spread. When X Is zero, a becomes C,H; thus, 

C, should be of the order of unity. . Og is calculated from o» as shown earlier, 

In the actual pool arm (X a 0) , the concentration thus becomes 
approximately 

Q/UC 3 C 4 H 2 g/m 3 

which assumes that there is homogeneity of concentrations to a height 

C-H and width C.H, This will be the highest concentration in the area 
3 4 

except in the immediate vicinity of the pool surface where the concentration 

3 
will be the saturation value of PM/RT g/m , 

In practice it will be desirable to express concentrations in 

3 9 3 633 

ng/m' by multiplying by 10 or ug/m using 10 or mg/m using 10 '. 

If parts per million by volume are desired, the concentration in 

,3 6 

g/m should be multiplied by 10 and divided by the substance vapour 

3 
density at 1 atm. pressure (g/m ) which is P M/1T, or at 25°C or -288K, 

or 42M where M is the molecular weight (g/mol) . 

It is likely that a is dependent on sampling time - which is here 

f 

.24 hours. For shorter sampling times the plume will be "sharper 1 ', o 

y 

smaller, and local concentrations higher. One could postulate a 
relationship of the form 

o ©< (sampling time) 
7 

where n may lie in the range 0.3 to 1.0. 
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2-4 Linking Evaporation and Dispersion 

For conceptual purposes, it is convenient to view the evaporation 
process as occurring as illustrated tear Her is Figure 2.1, in which there 
is a volume of well mixed air overlying the spill. This air is replaced 
at a rate determined by the wind speed U, If the evaporation rate is 
Q g/s, then the well mixed concentration above the spill will be 
Q/(UHh) g/m . The saturation concentration immediately at the liquid 
interface is PM/RT g/m . The blanketing effect is thus quantified by 
the ratio of these two concentrations, namely 

QRT/(UHhPM) 
But 

Q - K APxM/RT 

thus the ratio becomes in the case of a pure liquid 

K A/(TMh) 
e 

This is the dimensionle-ss ratio of two volumetric flow rates 
(m/s), the diffusion or evaporate rate from the interface, and the 

rate at which the vapour is swept from over the spill. Now A will 

2 
typically be the same order of magnitude as H ; thus, the ratio becomes 

K H/(Uh). K is typically O'Ol m/s, U is 2 m/s (7 S 3 km/h) ; thus, the 

e e 

ratio is unity when H/h is 200. The implication is that when the 
vapour passing over the spill is confined to a vertical height h which 
is l/200th of the spill diameter, then there will be blanketing. It is 
believed that vertical mixing is usually more rapid than this and that 
typically over a 10m spill the vapour has diffused vertically to of 
the order of lm; thus, H/h is 10 and the dimensionless ratio is 
0.05 or the evaporation rate is reduced by some 51 due to blanketing. 
Only in cases where wind speeds are very low, vertical mixing is slow. 



- 15 - 
and the vapour is negatively bouyant, will blanketing be severe. 

From the viewpoint of safety, the worst situation occurs when 
there is free evaporation and no blanketing; thus, in this work it is 
best to assume that the evaporation rate is not reduced by this effect, 
This will generally lead to a slight over estimation of concentration, 
It is unwise to assume that blanketing increases safety. 



- 16 - 



2,5 Vapour Pressures 

A substance's vapour pressure is a function of temperature alone. 
Values can be obtained from handbooks such as The Handbook of Chemistry 
and Physics or The Chemical Engineer & Handbook . Normally, the units 
reported are atm or mm Hg (torr), which should be converted to Pa 

1 atm m 101325 Pa 
1 mm Hg = 133 Pa 
If no vapour pressure is available, an estimate can be made from 
the substance's boiling point T_(K) using the Clapeyron relation and 
Trouton's Rule, namely 

P E = P B exp [10-6 (1 - T B /T E )] 
where F is the vapour pressure at environmental temperature T (K) and 
P g is the vapour pressure of 1 atm or 101325 Pa at the boiling point 
TgflC). Note that both temperatures are absolute (ie. °C + 273). 

For example, toluene boils at 110°C (383K) ; thus at 6°C (279K) its.a 
estimated vapour pressure is 

101325 exp [10-6(1 - 383/279)] 
= 1950 Pa or 14 -6 mm Hg. 
The actual value is 10mm Hg, 

This estimation procedure is approximate and tends to over- 
estimate vapour pressures and thus errs on the side of safety. It is 
applicable only to liquids, not to solids. 

If vapour pressure data P is available at a temperature T other 
than environmental conditions, then this equation can also be used, since 
if it is assumed that 

P D = P B exp [10-6 (1 - T B /T D )] 

then 

V P D * £XP [1°-6T B (1/T D - 1/T E )] 
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For example, the vapour pressure of toluene at 32 °C (305K) is 
4 Oram Hg; thus, at 6°C the ratio of vapour pressures is 

P„/P n - exp [10 -6 x 383 (1/305 - 1/279)] 
- 0-29 

thus j P is estimated to be • 29 K 40 or 11 "6 am Hg (actual value 
IQnim Hg)', 

It is preferable to use data as close as possible to environmental 
temperatures . 

As part of this work a literature review was undertaken of PCB 
vapour pressures, 'The results are given in Table 2.1, which contains 
estimates of vapour pressures of PCB' isomers and Aroclor mixtures. Full 
details of the estimation procedure are given in a separate publication 
(Mackay et al., 1982), 
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Table 2.1 
PCB Vapour Pressures and Boiling Points 



Aroclar 


Mol Wt. (a) 


BP Range °G 




Vap , Press . 


at 25°C^ 












atm* 


Pa 


1016 


258 


325 


- 326 


5.3 


X 


-7 
10 ' 


.059 


1221 


192 


275 


- 320 


8.8 


X 


■» /\™6 
10 


0,89 


1232 


221 


290 


- 325 


5.3 


X 


10~ 6 


0.54 


1242 


261 


325 


- 366 


5.3 


X 


7 (c) 
10 


0.054 


1248 


288 


340 


- 375 


6,5 


X 


7 (c) 
10 


0.020 


1254 


327 


365 


- 390 


1.0 


X 


10~ 7 


0.010 


12-0 


372 


385 


- 420 


5.3 


X 


io" 8 


0.005 



(a) 0. Hutzinger, S. Safe, V. Zitko, The Chemistry of PCB's (CRC press, 1974) 

(b) EPA Report 440/4-79-029 , Water Related Environmental Fate of 129 
Priority Pollutants 1979. 

(c) We believe that these values are suspect. Better values are probably 



1242 
1248 



8.0 x 10~ 7 atm - 0.08 Pa 
3.0 x 10 atm - 0.03 Pa 
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3, F1TTXMG TO ESBEEMENTM, DATA 

The model was fitted to the experimental data provided if the 
* Minis try in the report, "Poly chlorinated BIphenyls Survey at Smithville" . 

'The approach used was to test various values of the parameters C^, 
C-, C~ and C, until a satisfactory fit was obtained. Other forms of 
correlation were also tested but only the final results are given here. 

The best values were: 

CL - 0.15 C 2 - 5.0 C 3 - 0.05 C 4 - 1.0 

When fitting the data, the procedure used was to minimize the error 
expressed as the absolute value of the logarithm, of the ratio of 
calculated and experimental concentrations. 'This avoids giving undue 
weight to high concentration points. The complete set of calculated and 

experimental data are given in. Table 3.1 and illustrated as a function 

of distance from source X, windspeed U, temperature, T and angle, 6 versus log of 

concentration ratio in figures 3.1 to 3.4. These graphs show the dependence 

of the ratio of calculated and experimental concentrations as a function of 
each variable. If the correlation, was biased, for example, overestimating 

the effect of distance,, this would result in points lying off the horizontal - 

possibly at an angle. The figures show no such tendency for any variable; 

thus we conclude the correlation is satisfactory in its treatment of each 

variable. 

It was found that a satisfactory fit at long distances could only be 

3 

obtained if a, constant background concentration of 0*75 n.g/m 'was added 

to the calculated values. 

'The exact area of the '''spill'"' was not 'known, in this case; thus, 
an estimate was 'made that it was -,.240 B^ but not as a single, pool. It was 
regarded as consisting .of numerous pools of a total area of 24Q 1 m z but 
spread over an area of 7000 m2; thus only 4% of the actual ground area 
was PCB covered . 

It, is apparent from the Table and Figures that there is a high, 
inherent variability in the data. This, it is believed,, is due 
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primarily to errors in windspeed and direction, the wind data being 
measured some distance from the site and averaged over 24 hours. 

It was found that 901 of the model estimations were within a factor 
of 4*0 of the experimental values, which is regarded as fairly successful. 
It is doubtful if any model can substantially improve the fit of these 

The value of C, of 0.15 is double the "expected" value of 0.07 

obtained from dispersion charts. This may be due to ground induced 

turbulence. For example, there are buildings and trees in the area of 

height 5 to 10 m; thus it is likely that, at a distance of 100 m, a q 

z 

increase of only 7 m is too small, the actual indicated increase of 15 m 
being in reasonable accord with the dimensions of the "objects" and 
ground undulations. 

The value of C_ of 0,05 is equivalent to 1/20 of the pool diameter. 
It is important only in the immediate pool area, 

C, being a coefficient of unity on H can be ignored as a parameter. 

C ? is 5.0 m/s reflecting the "velocity" of increase of the perimeter 
value of a . Under Class C conditions (6 m/s) between 100 and 1000 m, 

y 

C Increases, according to dispersion charts from 8 to 110 m, i.e. by 102 m. 
This corresponds to a transit time of 150 s and thus to a horizontal 
velocity of 0.7 m/s. But this is for a shorter sampling time of typically 
of 0.5 hours, whereas the present time was 24 hours. 

We postulate that o depends on sampling time, the above data 

y 

suggesting a relationship of the form 
if &&> (sampling time) 

or BL ^5.0 (sampling time/24) * 5 

0.5 
~1,0 (sampling time) 
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The "spreading velocity" C for short sampling times is thus 1 m/s which 
.seems physically reasonable. 

It oust be emphasised that these values are only estimates which 
must be improved by further determination. The form of the equations 
also requires validation. In essence, the work is regarded as a step 
towards a validated procedure. 

As a final note, the procedure developed here obtained an estimate 
of the evaporation rate from the Smithville site of 2.0 x 10 g/s. 

A copy of the computer program used to generate these results is 
appended. (Appendix I) . 

A "calculation manual" was developed and is given in the Appendix. 
The "manual" was used illustratively for the Isabella Street School incident 
discussed in the next section. 
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Table 3,1. Experimental and Calculated Concentrations at Smithirille, Ont. 



»iea no. saapl« 

^ ^—duration (h). 



24.0 
24.0 

?fc.il 



<<4.0 
2*. 
74. 



2*.0 
24. 
24,0 



24.0 
24.0 
24.0 



24.0 
24. 
24.0 

aw . 9 



24.0 
24.0 
24.0 
3 fli . i7> 



disc from 
i. a 

1*0 

L*a 



9 



wind sp*«d 
flna/M 



t«mp«r*ttir# 
f°CT 



cemcs axpd 
fna/m 3 ) 



150.0 
55.0 



9.0 
10.0 



1 3.0 



284*0 
660.0 

50. Q 



1. 9 
1.0 
l.Q 

U4- 



100.0 
IHO.O 
120.0 

tTn.li 



17.0 

1 3 . 
10.0 

1 1,1 II I 



6 . 
4.0 

t ,0 

-' , n 



3.2 

3.7 
177.0 

?i .n 



i. o 

l . 
1.0 

1 . a 



165.0 

it o. o 

170.0 



15.0 
17.0 
16. 

i?.n 



3.0 
13.0 

7.0 

i i . n 



l> o 

1.0 

1. 

1. I,i 



30.0 
oS.O 

45. 

' 5 f - ■ a 



9.0 

8.0 

13.0 



1 0.0 

1 §■ 

ia.o 
1 i -" 



14.9 
410.0 

53.0 
102.0 



395.0 

1 160.0 

935.0 



1. Q 
I. 
1. f 



105. J 

70.0 

125.0 

too , a 



11.0 
11.0 

19. o 

1 2 . 



13.0 

11.0 

fc.O 

n.g 



1. Q 
1. a 



125.0 

lab.O 
40. 

i ~i e .. a 



15.0 
20.0 
12.0 
11.0 



13,0 

13.0 
7.0 



AVERAGE CACT = 



405. 12 



AVERAGE CCAL= 



1 193.0 

900.0 

29.0 

78 . 



522.0 

1 15.0 
533. O 



conca. calc. 
tng/«3) 



124.9 
271.5 

10^.4 



37.0 

78.3 

117.0 

77. A 



71.1 
189.1 

96. 1 

i 'ny.a 



234.5 
319.9 
331.6 

i •'h. a 



213. 4 

215.4 

98.7 

13 .6 



180.7 
147. 1 
170.4 
136.0 



24,0 
24. J 

24.0 



24.0 
24.0 



1. 
1. 

1. 



120.0 

35.0 

I I n. 



9.0 

10.0 
23.0 



s.o 

1 3.0 

9.0 

— a^Ji. 



1. 
1. J 

i . a 



3 0.0 
ICC.O 
105 .0 



13.0 
10.0 
10. 



4.0 

6.0 

3.0 
3.0 



89.0 
970.0 
195.0 
258.0 



74. 1 
247.0 

52.7 
1 1 3 .Q. 



140.7 
287.0 
129.4 

1 3 . n 



108.5 
166 .8 
100.9 

ao.o 



«£4.0 
24.0 
24. 



1. 

1 . 
L. 



1 1 0.0 
tfc>C .0 
10". 



2b. 
17.0 

ie.o 



7.0 

13.0 

7.0 



44.0 

OS. 

1S0.0 

2.26.0 



107.4 
154.3 
125.9 

1 =19. 1 



^4.0 
24,0 

24. J 



1 . 
1. 
1. 



of. .0 
155. 
1 1 5 . U 



9,0 
8...0 

2 0.0 



10.0 
I 5.0 
1 1 .0 

. i s — a - 



498.0 

184.0 

81 .3 

2 8 . 3 



216.0 
218.3 
1S2.3 



24. o 

24.0 
24. 



1, 
l.Q 
1.0 



1 (. . 
1*5.0 
13C.0 



11.0 
19.0 
12.0 



t t .0 
6.0 
9.0 

13.0 



23. I 

34.0 

270.0 

100.0, 



146.6 

90.6 
136. 1 



?*.0 
24.0 

24.0 



1.0 
1.0 

1 . o 



10b. 
50.0 
45. 



£0.0 
12.0 
11.0 



13.0 
7.0 

9.0 



502.0 
t 18.0 
321.0 



137. I 
165.7 
200.8 



AVERAGF CACT = 



J34.57 



AVERAGE CCAL' 



159.15 



i C . M 



u a«4- 



24. o 
24.0 
24.0 



1. 
1. 
1.0 



145.0 

0.0 

2C .0 

■ oc a 



10.0 

23.0 
17.0 



1 3.0 
9.0 
6. 



61.0 
48.0 
91. S 
4S . fi 



182.9 
183.0 
151.8 
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24.0 
24.0 

44,3 



1.0 
1.0 
1, 



aO .0 

10.0 

5.0 



10. 
10.0 
15.0 



6.0 
3.0 
3.0 



54.4 
20O.O 
I 19.0 

i fta.,i 



141.6 

133.7 
123.2 



> 



24. 

?4.0 
2*.0 



1,0 

I . o 

1.0 



ICO 

20.0 

170.0 



1 6 . 

12,0 

9 ,0 



7.0 
1 1 .0 
1 0.0 



129.0 
139.0 

31.5 
lAA.J, 



168. 1 
242.8 
140.1 

aa o.a. 



"4.0 

;4 . o 



l.Q 
1. 



75.0 

5.0 

95. 

W>0. J 



1 3a 

20.0 
11.0 
.1.0 



13.0 

1 1. a 



166.0 
169.0 

136 .0 
laa. o 



347.0 
220.6 
224.6 

IJ7.0 



24.0 
2».0 
T4.0 
^.. . 



t.C 
1. 
1.0 

-I .. <,? 



Cf .0 

2C0 

35.0 

I ft I,. . J 



19.0 
12,0 

15.0 

i2.a 



6.0 

9.0 

12.0 

7 . a 



187.0 

136.0 

1 73 7 . 

ia.A 



t 44.5 
207.5 
2<S3.3 

IQrt -7. 



24.0 



1. 



11.0 



9,0 



123.0 



187.8 



AVERAGE CACT= 



AVERAGE CCAL= 



24.0 
24. 



t. 
1.0 



70.0 
135.0 



9.0 
10.0 

7 1 \.Q 



e.o 

1 2.0 



136.0 

69.0 

5— Z_ 



ITT, 
19 0.0 



2 4. j 

24. J 



I C . 
100. 
160. 



17. 
13.0 
I O. 



o. 

4.0 

6.0 



26.4 
28.7 

IO.0 



114 
103 

100 



35. 3 
dO.O 

30.0 



15.0 
2«i. 

1 7.0 



3.0 

7.0 

13.0 



23.5 

42.4 

335.0 



99 
1 24 

359 



24 .0 
2 4. J 

:4.o 



1,0 

1. J 

1, 



nG. 
1 10. J 



12. 

•i.a 



i 1 . o 
i o.o 

1 £.0 



105.0 

lrt.3 

309.0 



219 
172 
360 



24. ) 
24.0 
24.0 

■ a 4, a 



i ... o 

1. G 
1,0 



75.0 

175.0 

t CO 



20.0 
t 1 .0 
1 1.0 

I -i.e. 



1 1 .0 
13.0 
11.0 
— 4_ 3- 



404.0 
253.0 
256.0 
iO.-Q 



185 
16 7 
249 



.74.0 
24.0 
2*. Q 
■ £» . .-> 



I. 
1. 
1. 



4S.0 
8E. 



t 2. 
IS .0 
2 0.0 
\ 3 . 



O.O 

12.0 

13.0 

— *^&- 



114.0 

272.0 

112.0 

— 1-4-, 



187 
257 
206 

-A3* 



14. J 



I. 



1 1 .0 



9.0 



90.2 



AVERAGE CV 
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Table 3.1, (continued) 



— — 

J6S.0 

TO i 
5fl,Q 



9 

i 



.2 4 « 

24,0 

24. J: 

-P4,0 



75, 

75.0 

7-..0 



12*0 
9,0 
85,0 

13,0 



11.0 

1 0,3 
I 5 , 0' 



5,8 
2.4 

10,4 



35.1 

17, 6 
37.5 
»3,0- 



24.0 
24* 

2 4,0 

3 a . a 



7 5.0 

7 5. 
75, 
7 E . 3 



1 2 0,0 

6 5.0 

l o „ a 



20.0 
I 1 . 0: 
I 1 . 

10,0 



I 1.0 

13.0 
1 1 . 



7.5 

39.8 
7.6 



38,6 
.22.6 
28,5 
44*6. 



24.0 

24,0 
24. 



7 5, 
75, 
75. 
75^-0 



5.0 
13.0 

30,0 
17EW 



12.0 

I 5, 
20'. 0' 
1 I ■ g 



9.0 

i s-« a 

13.0 



13,7 
19,9 
18.6 
-1-.-6- 



30,1 

43,3 
43, 1 

a. 7 



§ 24.0 

AVERAGE CACT= 



7b, 



90'. 
AVERAGE CCAM 



11.0 
28.09 



9.0 



12.2 



20,9 



2 4.0 
24.0 
34.-0- 



27 0,0' 
2 70. 



1S.Q 
170.0 



9 . 0' 

i a ,. o 



0*0 

13,0 



3,3 
2,0 



4.1 
3,2 



24,0 
.24,0 
24,0 
2 ft.. 



f7@.,o 

270. a 
27 0,0 
JAB fc, 



4b. 

45.0 

105. 

S^ i 



1 7 . 
13.0 

1 . 

34, a 



6.0 
4,0 
6.0 



1.0 

1.4 
1.5 



4.1 
3,4 
2.7 



24.0 

24,0 
24, 
3 ' i ., Q 



27 0.0 

27 0,0 
270, 
. 57Q.C 



25. 

35,0 

6,0 

-ia-&^o- 



17, 

1, 6 . 
12.0 
9,9- 



13,0 

7,0 

11.0 



3.0 

0.5 

3,3 



7.4 
4,6 



24.0 

,2 4 . 
24.0 
-4-4. 



24.0 

24.0 
24.0 



27 CO 

.2 7 0, 

27 CO 
T70, o 



270. 
27 0,0 
2 70. 
■WU*. ■ ■ 



7C0 
20.0 

ise.o 

- 4&»-9- 



2 0, 
11*0 



I 0.0 
5*0 

I 0. 



19,0 
1 2 . 

15,0 



15,0 
1 1 .0 
13.0 

11 . 



6.0 
9.0 

e 



ii: 



5.e 

3,2 
7.9 



5.0 

2,3 
1,4 



5.7 
7.0 
3,6 



5,0 
S.O 
7,4 



2 4,0 
2*.Q 



2 7 0. 
270, 



1,75,0 
90. 



.12.0 
11.0 



7.0 

9,0 



0,6 
2.1 



1,8 
3,5 



.ncrmrg r * /-t . 



4-W--EAr,£ CCAV?_ 



2*.0 

2 a , -a 



35 0', 

as s. o - 



28,0 
37„9- 



9, 

19.. a 



.fi.O 

13.0 



2.0 

0,4 



2,8 



2 4.0 
2 4.0 
2 4 , 



3 50. Q 
350. Q 
35 0.0 



I S. 
2,0 

a . o 

A3, SL- 



23.0 
17,0' 
13,0 



9,0 
6.0 
4.0 



1.3 



4,5 
3.4 
2.7 



24.0 

2 4,0 

24.0 



35 0.0 
350, 
350. 



S.O 
1,3,0 
18.0 



10,0 

15,0 
2 6.0 



3. D 

3,0 

7,0 

7,0 
11.0 

10,0 



3,6 
1.3 
0.2 



1:1 



2*. 
24.0 
2*.0 



35Q, 
350. 
250 , 



6.0 

ja, 

1S2#0 



16,0 
12.0 

9,0 



. 1 

0-1 
0*4 



3.5 
3,7 
,2,2 



9 



2 4.0 

24.0 
24.0 



2 4. J 

24,0 
24.0 



350.0 
35 0.0 

3 50, 
-M4 



9.3.0 

23.0 

77.0 

-i98.0- 



350, C 
350, 
350. 

iaa« o- 



53, 
3B, 
53,0' 



11.0 

- * t--=o. 

1 9 . O 

1,2.0 
1S.Q 
2 0,0 



1 8,0 

1 1,0 
13.0 



9.-0 
13,0 
13,0 



2.6 
1.0 
4,5 



1,7 
0.6 
0.9 



3.9 
4*8 

3.4 

a. a 



2,7 
3.3 

4*1 
5.7 



24,3 
24.0 



350. 
350.0 



I 42 , 
47.0 



A wsPAgg cefti..- 



1.2.0 

11.0 



7.0 
9.0 



0.8 
1 • 



1 .5 

3.1 



24.0 



5 70- 

-7n-.-o- 



i os. a 

— 5Q..0 



9.0 

i a . a 



o! 



i^_ 



1 .4 

2.0 



24,0 
24.0 
2 4.0 



57 0.0 

570, 
57 0.0 

axe o 



O5.0 

75,0 
75,0 

i s i a 



23.0 
17,0 
13.0 



9.0 
6,0 
4. 



0,4 
0,5 
0.6 
1 «f 



1.0 
1 ,3 
1.3 

i - 1 



2 4.0 

.2 4 . 
24,0 



570, 
57 0. 
5 7 0. 



14 5,0 

as, a 

115.0 



1 7. 
Id. 
12.0 



12,0 

7.0 
11 ,0 



0.3 

„ 4 
0,2 



0.9 

1.3 
1,3 



24.0 
2*. 
24.0 



vt c „ o 

570, 
5 7 0, 



170.0 

17 0,0 
100.0 

_n,-a--_ 



a.o 

13,0 

20. O 



1 5. 

1 S.O 

I 1.0 



US 
0,<£ 

o.a 



i ,a 
1.2 
1 .1 



24.0 
24.0 

24.0 



570. 
57 0,0' 
57 0.0 
-ST.a^O- 



175. 0' 
130,0 
115.0 

1 1 a . a 



11.0 
19.0 
12.0 



11.0 
<6. 
9.0 

H3U43 



0.4 

0.6 
0,6 



1 .2 
0,9 
1.3 
1 . 1 



24.0 

2 4.0 

24,0 



37 0,0 
570,0 
37 0- 



OO.O 
65.0 
30.0 



20. 
12.0 
11.0 



13,0 
7,0 
9,0 



0.4 
0.8 
0.9 



1.4 

1. . 5 
1 .8 



AVERAGE CACT= 



AVERAGE CCAL= 



c 

G 

G 



34-, 0- 



24, 
2 4, 
24.0 



2*0. 
2 4C. O 
24 0.0 



30.0 
17 5.0 

155,0 



10.0 

23.0 

,17,0 



13.0 
9,0 
6.0 



12.0 

0,7 
0.5 



7.0 

O .8 



,2 4 ■ a 

24,0' 
2 4, 



2 40, 
24 0.0' 
24 0, 



9 5.0 
1 65. 

1,70.0 



1, 0.0 
10.0 
15.0 



6.0 
3*0 
3.0 



4,0 

0*8 

0.9 



3,4 
2.1. 
1 ,2 



24.0 
24. 

2*. 3 



2*0, 
2*0, 
2*0, 



165.0 
16 5. 

5.0 



lte.Q 

12,0 
9. 



7,0 
11.0 

1 0. 0' 



0.5 

o . a 

15.4 



1 .3 

2 ,5 
5.6 



2 4. 
2 4.0 
24. 



240. "J 
24 0', 
240.0 



1 1 0.0 

ItsO.Q 

eo.o 



13.0 
20.0 

11,0 



1 fi • 
11.0 
13.0 



5,7 
1 .0 
5.7 
4.1 



2 4.0 
24. 

24.0 



,24 0. O 
2*0. 

24 0. 0' 



i be . o 

165,0 
150.0 



19.0 

12,0 

15, 



6, O 
9,0 

13.3 



0.2 
5.3 

0.4 



1 .1 

2,3 



?4, a 

24.J 



240*0 

2*0. 0' 



1 5. 
110.0 



12,0 
11,0 



7,0 
9.0 



14.3 

9,0 



5,2 
3,6 



AV£r*AS&--CAC-T- 



- AV6fM' , ;£-^C CA 4H = 



= 23 - 



table 3.1. (continued) 



I 10 
10 


2*.3 


1 000. 

1 .1:C|:(. n .V 


- 135.3 

- - "lil 


9.0 
1 ft t « 


e.o 


1.6 

i .n 


0.9 


= — 


10 
19 
10 

Ma . -. 


24. v* 
24. 

24.0 
34. J 


1 30 0.0 

1 ooc.c 

1 300. C 

i o * e ■ it 


125.0 
105. e 
1 05, a 


2 3a 
17,0 
13.0 

1_«, , p 


9.0 
fe.0 
4 .0 


0.6 

0.7 

o.a 


0.3 
0*3 
9.9 




10 
10 
10 

10 ■■ -■ 


24.0 

24. J 

24.0 

a i. . .i 


i ops, g 

1 00 0. 
I 000.0 


125.3 

175. 3 

1 IK. 3 

I 4 »- -? 


2£.0 
! 7.0 
IbtC 


7.0 

13.0 

7.0 


0,5 
0,2 

0.4 


9.3 
0.3 

O.a 




10 

10 
10 

_ui 


24. 3 

24. J 

24.0 

?»l 1 


i aoe.o 

1 OPS. o 
1 000. 

ijm^lv 


*5» tS 
1 40.3 
t € 0. 


V. 

8.9 

13.0 


lg.0 
I 5 . 


1 .6 
2.5 

t .5 


Ul 
1,1 
0.9 




10 

to 

10 
10 


24. a 
24. C 
24.0 


1 000.0 
1 30 C. 

1 OOC.Q 

i a a <■ ■ o 


30 .4 
1 *£. 

1 1- . 3 


11.0 
11.0 
19.0 


13.0 

1 1.0 

c. 


2.7 

0.8 
0*6 


1.2 

0,9 
0.3 




10 
10 
10 
I n 


24.0 
24.0 
24.0 
"a 1 1\ 


l out:, o 
l ooo. a 

1 90 0. 
1 ahn ft 


luii.O 

120.0 
35.0 

Mi ii 


1 5. 
29.0 

12.3 
111 


1 .1.0 

13.0 

7.0 


t .6 
1.2 
l.T 


0.3 
0.8 
1.1 




AVERAGE 


CACT» 


1. It? 


AVERAGE CCA.L* 


0.92 










11 
11 
11 


2 4.0 
2 4,0 
24.0 


1250.0 
1 250.0 
1 250. 

i a* o . a 


160.0 

5,0 
150.0 

i. j ?n . a 


9.0 
10.0 
23.0 


8.0 

13.0 

9.0 


9.6 
1.5 
0.5 


0.9 
1.1 

O.S 




1 t 
1 1 
1 1 
1 1 


24. 

24.0 
2 4.0 


1 250.0 

1 250. 
1 250. 

i 1 ■ani-i a 


130. 

70.0 
150.0 
ni f , fl — 


13.0 
10.0 
2u.O 
1 T . 1 


4.0 
6.0 
7.0 


0.8 

1 .4 

0.2 


0*3 
0.9 
0.3 




1 I 
I 1 
11 
1 1 


24.0 

24.0 
24. 0' 
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Figure 3. la. 



Concentration ratio as a function of distance from source 
for Sites 1-4, 
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Figure 3.1b, 



Concentration ratio as a function of distance from source 
for Sites 5-14. 
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Figure 3.2a. Concentration ratio as a function of windspeed 
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Figure 3.2b. Concentration ratio as a function of windspeed 
for Sites 5-14. 
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Figure 3,3a. Concentration ratio as a function of temperature 
for Sites 1-4. 
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Figure 3.3b. Concentration ratio as a function of temperature 
for Sites 5-14. 
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for Sites 5-14. 



200.00 



U) 
K5 



- 33 - 



4, AFP LIGATION TO ISABELLA STREET SCHOOL 

On October 3, 1980 a transformer at the Isabella Street School 

Thunder Bay, Ontario raptured and spilled approximately 130 litres of fluid 
containing approximately 70% FOB compounds. The transformer was situated 

near the east wall of the school at ground level, The spHI covered 

2 2 

approximately 100 m adjacent to the school and 10 m is patches up to 

25 m from the school, "The Ministry of the Environment monitored the con- 
centrations in the school yard .and building for several days after the spill. 
Using prevailing temperature and wind speed; and spill dimensions 

and sampling time supplied by the Ministry, the calculation procedure 
described in. ."Appendix I was used to calculate concentrations at various 
disptances from the school, The results are given in Table 4,1, 

Concentrations .measured indoors at the school are not calculable 
using this model. A fair estimate would be to assume that indoor and 
outdoor concentrations were equal, especially if there is rapid air exchange, 
In the event that the liquid was spilled in part indoors, an entirely different 
situation would apply. 

Table 4,1 Concentrations at Isabella St, School at various distances 

from the spill 



Conditions 

- 0,0° 

U « 4,2 m/s 
(15 Wh) 

T - 280 K 
f ^ 0,0048 
A - 100 m 2 



Distance from 

spill edge (m) 


Sampling 
Time (h) 


Concentration (ng/m ) 


at height m 


at height 1 m 





1 

i 


4200 


"560 


3 


2300 


1500 


10 


1 


830 


730 


. 20 


1 


402 


370 


23 


24 


140 


140 


100 


1 


39 


39 
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3 
The results show that concentrations of 1500 to 2300 ng/m are predicted 

at I distance 3 m from the spill edge at heights of and 1 m. These values 

are higher than those obtained experimentally by about a. factor of two. At 

3 
10 m distance , the concentrations of 730 to 830 ng/m are close to those 

obtained experimentally, By 20 m the concentrations have fallen below 

3 
500 ng/m for short sampling times (1 hour) - which is the worst case. The 

long (24 hour) sampling time values are a factor of approximately 3 lower. 

3 
By 100 m distance, levels of 39 ng/m are expected, well below any level 

of concern. 

It must be emphasised that the air flow characteristics in this region 

may be quite abnormal due to channelling between buildings, and flow over and 

around buildings. It appears that these flow "abnormalities" cause an increase 

in a in the region of the spill, i.e. vertical mixing is quite rapid. It 

z 

is probably impossible to quantify such flows accurately £ thus the only approach 
appears to be to use the calculation procedure as a pessimistic estimate and 
interpret the results in a toxic sense with an awareness that the concentrations 

may be overestimates by factors of up to 5. There may be periodic air gusts 

3 
with very high concentrations even as large as 13,000 ng/m (the value over 

the pool) but with others containing negligible PCS, 

In general, this test of the equation is taken as being satisfactory 

in that the calculated concentrations are within factors of 2 of those 

observed. It is puzzling, however, that some very much higher concentrations 

were not obtained close to the transformer. 



j 
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-5. DISCUSSION MB CONCLUSIONS 

A calculation procedure has been developed and illustrated in which 

air concentrations can be calculated in the vicinity of spills of evaporating 
liquids. It is emphasised that the contribution here is to bring together 
the component processes into an overall package which it is believed will 
be useful to the on-scene 'response team. There are several areas of 
improvement possible as follows. 

1, The vapour pressure estimation from boiling point applies only for 
non-polar organics. A listing of vapour pressures at various temperatures 
for "candidate''" substances would be valuable, 

2, The evaporation mass transfer coefficient depends on pool size and on 
the evaporating substance, A better correlation is possible. 

3 , The values of a and a depend on atmospheric state and should thus 

be taken as "typical values"'. Under stable conditions when the plume is more 

confined the procedure will underestimate concentrations at ground level on 

the centreline. Some judgement is thus necessary when selecting values under 

stable conditions. For example, at night under a clear sky,, a and a_ 

y z 

may be lower by a factor of 2 or 3. 

4, The pool temperature may be higher or lower than ambient depending on 
radiation and evaporative cooling. Again, judgement is necessary, 

5, This analysis does not apply to plumes which are sufficiently concentrated 
that they have positive or negative buoyancy. Of particular concerns are dense 
plumes (i.e. cold or with high concentrations of high molecular weight gas) 
which will spread over the ground at high. concentration. 
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6. There is no doubt that sampling time enters the calculation but the 
equation presented has little justification, A better equation should 
be devised, 

7. The vertical coef f icient a i is influenced by terrain "roughness", 

probably increasing by approximately the magnitude of the size of the object. 

Again, judgement is necessary. There is a danger, however, that overestimating 

it will lead to underestimation of concentration - which is unsafe. The 

safest procedure is to underestimate a . 

z 

In summary, it is believed that the procedure developed has sufficient 
rigour that the results are credible, but the concentrations must be regarded 
as estimates. Several improvements are possible and could be explored. It 
is recommended that the procedure be tested in the field to gain an appreciation 
of its applicability. Is it too complex, or insufficiently complex? Such 
experience would be valuable as a guide for future work. The procedure 
also highlights the monitoring data which could be taken at spill sites 
which could be used to "validate" the procedure, Such data (and those 
discussed here) are an invaluable source of information about spill evaporation 
and dispersion and it is hoped that they can be exploited, as has been 
attempted here, in the interests of providing better response advice 
following accidental spills. 
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CCAL( 1 )*Q( 1 )*0( I )»E( I 1*M<I )*C8ACK 
CR< I )»CCAL(I )/CACT( I ) 
CRL< I)=ALOG(CR< 1) ) 
- fcOMU I > s ASS ( G3 L( 1J, 1. 



CR( I )*A8S(CCAL II J/CACTU J ) 

Ci»AT( I )=VP( I l«WM»lE09/(R*TK(t) )*-CBACK 

CQC I )»Q( 1 )/(W5MSCI )*0IA*01A»C3>+CDACK 



10 CUNTlNUE 

nn ia .i-ai.iA 



TCACT=0 
TCRL=0 .0 
TCK1.A*0.0 



DU 30 lal.312 
IF(NN( I ) .EG. j I G3TQ 40 
GUT J 2 9 
/»C CilNTtMUli 



TCACT= TCACT*CACT ( I ) 
TCCAL»TCCAL<-CCAL( I ) 

TCPL=TCHL+CKL( I > 

1-Tf -r\.f 4C4k*<.I ) 



WdD( JL ) = W0< JL) 
XXt JC)=X ( I } 
*«bKH< JU)=WSKH{ I ) 
TKKt J L l -TK. ( i ) 



.Uu( JL )=WD< I ) 
SSntPt JL) = 5JIK ( I ) 
T"LU t JL)=Cf-L( I ) 

■I J- J Jit 1 



JL=JL+ 1 

»^ JTC{ ft. 102)NM< 1 ). SO ( I > .X< I 1 .*DI t ) ,wSKH( I ) ,TCt 1 J.CACTU J .CCALt I) 

■t.-.i.it imi^- 



30 CCNT INUF 
JK=J J- 1 

AVCACT«TCACT/JK 
Atiirii. Ah~TCC Aim,/ JJt, 



AVCRL = TCHL/ JK 

AVC^LA=TCRLA/JK. 

WW 1TEC o. 103 ) AVC ACT, A VCCAL 

J J= i 
20 CUNT 1HUF 
SCR = SCf»/'N 
I.DTT6 ( I, , in'iil S^T ■- - 



j 05 FORMAT ( '0' 
PETURN 
CND 



AVERAGE CR=».F1Q.6) 



- 38 - 



SAMPLE CALCULATION 



CA1XH1.AT ION PROCEDURE 

1 . Environmental condltloni 

H - *> m/a p ka/hr - 0.27B ■/■) 

°C T. - t. * 213 - %S^ R 



Hind Bj.ee d U 



Temperature t„ - I 

K - 0,0025 ll "" - O'OftTC ./ s 

i — - ~ 

2, Substance Properties 

Molecular velght H » Sj 2. g/mol 

Hole or volume fraction ■ - Oh; Co traction 

Density I - J ' sV, ■/*» 

Vapour preaaute P - • O'd'M a*_ 1?" (■»• below) 

Bote V apour Pressure <VF) 

la eonrert at* to Pa multiply by 101323 

1 ■ m ig to Fa "" *" 131 

If IP P. Is available at another temperature t °C then estimate the 
environmental Up P at temperature ^ C tram 

fy - r D expI10.6(C a + 2I9)(l/(t p * I7J) - l/fkg + *W»1 

wliere t R la the IioUtag point °C. 

If no 'UP la available but tin balling point t "t: is Inmni then 

P e » 101125«p(10.6(l - fig * J?l)/(tj * M3)) 



3. Spill Properties 

Vol uae spilled « - Q I H3, « J Hotel 1 ■? - ».J ft 3 

- J20 Igall 

- 264 USgall 

- 1000 Utnn 

- 6.1 barrel a 



Area of spill a - _ __ I Q Q m 2 

Itepth #f spill - ¥/A - _ fc> . <a I 3 ■ 

Disaster of pool II - I O ■ 

4. Estimate Evaporation Kate gj 

Evaporation rate q ■ I^APjift g/s 



ultere i ■ 1,3M 

T„ - t„ + 273 



Spill lifetime i ■ 10 fi Vd/Q a - S m ¥7jM P' ' a 

Note 11 spill depth D Is known 

I - 10 6 lhlRT/<K Pun) a - 



S. Estimate Concentratlona In. Area ot Pool 

|IJ HsalmuB) concentration at spill surface 

""I 

c: a - PH/If g/- 3 ' - Co •Ci'X g£g E f m 3 

(111 Concent rial Ion above pool 
C - 11/(0. 1)5111 Z ) 






*. lettagte DlBBsrslon Cpstflclsnti 

Desired distance X ■ 3 ■ 

ismpltal time: * - j hoars 

Vertical coefficient o - 0-0511 * 0.15s 



motliontal coefficient a 



III + S""s/U 



-/Q'7 - 



Angular coefficient o^ - 1B0 bM - T^PHlH deirtea 

(If o fl excaeda 180° use HJ q and recalculate 



O from 
ff 



bK 



gj^ 



7. Calculate Centreline Crowd Level Concentration at I 



Q/(sUo ]j o y ) 



5-3^70 g/a^ aSob ng/m* 

C alculate CoMcentratlon Of f -Ground and Of f -Centreline 

- 

Height I - I ■ 

Q 



Angle between Mind and subject I 



light factor r - exp(-s Z /(2o .*) 



degrees 



o - £T7i" 



•Angle facto* F g - eap(-o 2 /(to a 2 l * «»ptS»-»0> 3 |ti»« t | 
-_ /■ tJkk degress 



L<.' i; t.ei,.. : ■«■: ;:.;- ;, 



ir &t* 



j Sjfl O Jlg/m 3 



SAMPLE CALCULATION 



APPENDIX I 



CM.Rtti.AT (OH RKWBTOKi 

I, Environmental conditions 



Ulnd speed M 

Temperature l 
K - 0.002S M 6 



1 - 



■7b (i Win: - 0..2J1 mha 



,1, Substance Properties 

Molecular weight M ■ 3^ "3* !#■** 

Mole or vol one fraction ■ - Q*_Cp traction 

Density d - I j % g/c« 

Vapour pressure IP - • OOH? Pa |mms below) 

dote Vapour Preaaiir ja (VP) 

To convert situ to Pa multiply by 101125 

" an III W M " " 113 

If ¥P P la available at another temperature t ' C then estimate the 
environmental VP P g at temperature t^'C twos 

P • F exii|t0.6(t a + Z73){l/(t D * 273) - l/t« E * i»il'J 

where t la the boiling point °C. 

If no. VP It available but the boiling point t "c le known then 

t % •& 101]25exp(10.6<l - (!„ * 17})/ (t,, + 213)) 



1. Spill Properties 

Volume eplllei V - A . Q% a? Hotel 1 m? - 35.3 © J 

■ 220 Igali 

■ BH USiall 

- 1000 lltraa 

- 6.3 barrels 

Area of spill A - f Q t> m 

Pepth of spill - V/A - fa t QD t 3 ■ 

Diameter of pool I - I O ■ 

4 , gafjawte Evap oration Rate 

Evaporation rate q ■ KJkBM g/s where I«1.3H 

Kt t T E - t |E * 273 

- n-isy/o pi. 

Spill lifeline % m 10 6 Vd/Q ■ - 5 " %7rWP' ' j 

Mote II iep.fl 1 depth D ta known 

1 - 10*lklkT/|K-PxH) a - " a* 



5. Estimate Concentratl_qna_ln Area g£ f rot 

(I) Hail nun concent rat ton at spill surface 

(It) Concent rat ion above pool 



!-yi/o E/B 3 



i. Eatlmate Dlapatalon Cwff iet>„t« 

Daaired distance I ■ 3 

Sampling tiM S - 



I 



Vertical coefficient a) - O.OSH * 0. 1SI - . Q *°l S* u 



Hotlxontal coefficient O - H + S *X/U 



ma 



Angular coefficient o § - WO o /(aX) - JLOf'*/ degre 



(If o exceeds 1B0° use 180° and recalculate 



"8 
f 

P 



^JS . 



J. Calculate Centreline Ground Level Concentration at 11 

c c - Q/irttyy . 3-3 X/o" g/J = a3op 

■• Calculate- Concentration Off-Ground and Off-Centreline 

Height Z - I. ■ 

Angle between wind and aubjoct S - Q <J, 



ng./m 3 



Height factor F 

1 

Angle factor P. 



mrirrfOto*! 



Q • S'7S" 



«ipl-* 2 /(2o e ! ) * eap(-<B-i&n) 2 /(2o Q 2 > 



/./IS" 



degreea 



Concentration C 



/■if <**» 



. & 



■ SO O .pg/m- 






Hi in i ir ii ii mi in 

* c ltT3bDDDD l DTS7h* 



